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ABSTRACT 

Context. Diagnostics of MHD waves in the solar atmosphere is a topic which often encounters problems of interpretation, due partly 
to the high complexity of the solar atmospheric medium. Forward modeling can significantly guide interpretation, bridging the gap 
between numerical simulations and observations, and increasing the reliability of mode identification for application of MHD seis- 
mology. 

Aims. Determine the characteristics of the fast MHD sausage mode in the corona on the modulation of observable quantities such as 
line intensity and spectral line broadening. Eff'ects of line-of-sight angle, and spatial, temporal and spectral resolutions are considered. 
Methods. We take a cylindrical tube simulating a loop in a low-/? coronal environment with an optically thin background, and let it 
oscillate with the fast sausage mode. A parametric study is performed. 

Results. Longitudinal structuring of the intensity modulation is obtained, set by the nodal structure of the radial velocity. The mod- 
ulation is strongly dependent on the contribution function of the spectral line. Under the assumption of equilibrium ionisation, the 
intensity variation can be very low (< 4 % for Feix 171) or significant (35 % for Fexii 193). Most of this variation disappears when 
considering the radiative relaxation times of the ions, due to the fast timescales of the sausage mode in the corona. Regardless of 
the ionisation state of the plasma, the variation of spectral line broadening can be significant, even for low intensity modulation. The 
nature of this broadening is not thermal but is mostly turbulent. This places spectrometers in clear advantage over imaging instruments 
for the detection of the sausage mode. The modulation of all quantities can considerably decrease with the line-of-sight angle with 
respect to the perpendicular to the tube axis. The spatial and temporal resolution are the main factors aff'ecting modulation, erasing 
longitudinal structuring when these are on the order of the mode's wavelength or the mode's period, placing high constraints on 
instrumentation. Significant variability in all quantities can still be obtained when viewing at an angle of up to 30°, with pixel size 
resolutions up to 1/3 of the mode's wavelength, or temporal resolution of one fifth of the mode's period. Modulation is only weakly 
dependent on spectral resolution due to the mode's inherent symmetry. 
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1. Introduction 

During the last decade solar physics has witnessed a boom of 
observational reports on wave-like behaviour in all layers of 



the solar atmosphere (see the recent review by |De Moortel & 
|Nakariakov|20 1 2| ) . This is due, on the one hand, to the advances 
in imaging and spectroscopic instruments, allowing higher spa- 
tial, temporal and spectral resolutions, which have allowed the 
proper visualization of MHD waves with a rather large range of 
properties in terms of amplitudes, periods and speeds. On the 
other hand, theoretical progress points to the strong potential of 
waves for explaining the heating demands of the chromosphere 
and the corona in the observed spatial and time scales, as well 
as the important role they can play for the indirect estimation 
of local phy sical par ameters of the plasma (MHD seismology, 
Roberts et al . 1984; Uchida|p^7Ql |Aschwanden||2004t [Taroyan 



& Erdelyl|2009| ). This of course applies to all physical param- 



eters dic tating the behavi our of these waves: the magnetic field 
strength ( Nakariakov et al . 1999; Pascoe et al. 2007b ), th e mag- 
netic field geometry ( Verth & Erdely 1,^200 8 ^ ,Ruderman et al.| 
2008[[Wn Doorsselaere et al.|2009|), the density inhomogeneity 



(|Aschwanden et al.||2003| [Goossens et al.||2002| [20081 |Arregui| 



et al. 2007 1, the temperature structure, the response of the plasma 
to ener gy injection (poly tropic index, |Van Doorsselaere et"aL 
201 lb| ), the plasma JS parameter ( Van Doorsselaere et al.|201 la ), 



and the transport coeflftcients ( Nakariakov & Verwichte|2005a| ). 
Therefore, much emphasis has been put on the search for all 
kinds of MHD waves, each contributing in a more or less spe- 
cific way to the measurement of a specific quantity. 

The fast MHD sausage mode (m = mode), which we 
consider in this study, is compressible, characterised by axi- 
symmetric radial displacements of the plasma, having magnetic 
pressure and gas pressure acting in phase as restoring forces. 
Thus, this mode does not perturb the main axis of loops but 
its cross-section, producing intensity variations and spectral line 
broadening. 

The MHD fast sausage mode has been first de scribe d analyti- 
cally b y|Rosenberg|([1970),|Zaitsev & Step anov ( 1 975| ), [Roberts 
[et al.[ ( [1984| ), and [Cally[ |l986| ). Most observations of oscilla- 
tory phenomena interpreted as sausage modes have been made 
in events associated to solar flares, observed with radiometers in 
the microwave range. These have been most frequently reported 
using the Nobeyama Radioheliograph (NoRH, Nakaj ima et al. 
[I994j), with a cadence of 0. 1 s and a spatial resolution of 5'' - 10'' 
observing in microwave frequencies of 17 GHz and 34 GHz 



CNakariakov et al.'2003[ [Melnikov et al.|2005|[Reznikova et al 
2007 ; Inglis et al. 20081 The RHESSI spectrometer ( |Lin et al. 



2Q02| ), with a time resolution of 4 s and a spatial resolution of 
2.3'' in the 3 - 100 keV range, has also been able to observe 
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modulat ion characteristic of sausage modes in the hard X-rays 
range ( Zimovets & Struminsk y||2010| ). In observations with the 
LYRA radiometer aboard PR0BA2, which includes both chro- 
mospheric and coronal channels, Van Doors selaere et al.| ( |2Ql la| ) 
report two modes of oscillation in a single flaring loop, which 
were interpreted as slow and fast sausage modes. Few reports 
also exist in EUV lines with the EIS spectrometer of Hinode 
( Kitagawa et al. 2010|), w ith the Atmospheric Imaging Assembly 
of SD Q ( |Su et al.|2012| ) and with SECIS ( [Williams et al.|2Q021 
200 Similarly, few reports also exist in cooler (chromospheric 



and photospheric) lines, with ARIES ( [Srivastava et al.||2008 



with the Solar Optical Telescope of Hin ode ([Fujimura & Tsuneta 
| 2009), and with the ROSA inst rument ( [Morton et al.|2011t|Jess 
let al.|2012l|Morton et al.|2QT2l ). 

Most observations of fast sausage modes are linked to 
Quasi-Periodic Pulsations (QPPs) commonly observed during 
solar flares in microwave and sof t and hard X-ray ranges (see 
INakari akov & Verwichte[ [2005bl [De Moortel & Nakariakov 
2012, for reviews on the subject). As shown by Mossessian & 



Fleishman (2012), not only the sausage mode but also the kink 



and the torsional Alfven mode can lead to modulation of the 
source parameters (e.g. energetic electron distribution, magnetic 
field strength and column thickness along the line-of- sight) re- 
sulting in gyrosynchrotron radiation. 

Nakariakov et al.[ ( [2003] ) have shown that the period of 
the sausage mode in the corona is mostly determined by the 
length of the hosting structure and the external Alfven speed, 
thus leading to very short periods. They further show that non- 
leaky sausage modes occur in sufficiently thick and dense loops, 
thus greatly constraining their existence in the solar atmosphere 
([Nakariakov et al. 2012 ). It is therefore not surprising that most 
of the reported sausage modes in the literature come from ob- 
servations with high cadence instruments, mostly radiometers, 
and are mostly associated with flaring loops, where densities 
are high. However, by performing 2D MHD simulations, [Pascoe 



et al. ( 2007b ) have investigated the existence of sausage modes 



in long coronal loops with weak density contrast by extending 
the wave numbers into the leaky regime. In flaring conditions, 
they found sausage modes with periods from 5 s to 60 s of de- 
tectable quality in loops of length up to 60 Mm, thus supporting 
the interpretation from observations of QPPs. In this regime the 
period of the sausage mode is determined by the length of the 
loop and the external Alfven speed, therefore providing a seis- 
mological way to measure the latter and the external magnetic 
field. Very few reports exist, however, of sausage modes in long 
coronal (non-flaring) loops, a fact that may be linked either to 
MHD radiation (leaky modes) leading to a poor quality factor 
for such modes, or to the lack of adequate temporal and spatial 
resolutions. 



In a 3D MHD model [De Moortel & Pascoe[ ( |2QT2l ) ana- 
lyze the eff'ects of line-of- sight integration on energy budget es- 
timates in multi- stranded coronal loops subject to propagating 
kink modes. The modes strongly couple to azimuthal Alfven 
waves, to which all energy is eventually transferred, resulting 
in kinetic energy estimates from Doppler shifts 1 to 2 orders of 
magnitude lower than the true wave energy. This suggests that 
most of the energy could be observed in spectral line broaden- 
ing. 



Cooper et al.[ ( |2003a|b[ ) investigate analytically with a 2D 



model the modulation of intensity from propagating kink and 
sausage modes based on variation of the angle between the line- 
of-sight and the emitting structure. They further consider the ef- 
fect from curvature, and apply the results to explain propagat- 



ing intensity perturbations along a loop observed with SECIS 
(•Williams et al.i2001','2007). 

Gruszecki et al . (2012) have performed 3D numerical sim- 
ulations of (trapped and leaky) sausage modes oscillating in 
a cylinder, and study the resulting modulation of the intensity 
along a ray perpendicular to the axis of the cylinder. They con- 
sider the variation of the pixel size, thus simulating diff'erent spa- 
tial resolution during observations, and find that intensity vari- 
ation essentially disappears for resolutions on the order of the 
mode's wavelength. 

These previous works calculate the intensity by integrat- 
ing the density squared along the line-of- sight, and therefore 
ignore the radiative response of the plasma and eff'ects from 
instrumental resolution. Based on numerical simulations, syn- 



thetic observations of TRACE channels and Hinode/EIS by De 
Moortel & Bradshaw] p008|) illustrate the 'ill-posed' intensity 
inversion problem, which can even alter the period of oscilla- 
tio ns. Similarly, synthetic observ ations of AIA channels of SDO 
by Martinez-Sykora et al. ( 2011| ) show that significant contribu- 
tion to the emerging intensity can come from non-dominant ions 
of a specific channel. These results imply that not only density 
but also ionisation balance and emissivity need to be taken into 
account for forward modelling intensities. 

The application of MHD seismology strongly relies on the 
capacity to pinpoint the nature of the wave in observations. 
However, as analytical, numerical and observational studies 
show, mode identification is far from a trivial matter. In this work 
we choose to use a simple approach to characterise the observ- 
able properties of the MHD fast sausage mode in the corona. 
This approach can be used as a platform for constructing more 
complex models or can help understand more complex numeri- 
cal results, such as those coming from 3D numerical simulations. 
We assume that we have a cylindrical structure in a low-/3 coro- 
nal environment, which off'ers significant density contrast with 
the ambient plasma, and hosts the fast sausage mode. By syn- 
thesising the emerging intensity in specific coronal lines using 
the CHIANTI atomic database ( [Dere et al.[|1997[ 12009] ) we in- 
vestigate the eff'ects from varying line-of-sight angles, as well as 
eff'ects from spatial, temporal and wavelength resolutions on the 
modulation of observable quantities. By doing so we are able to 
predict the detectability of this mode with current instruments, 
as well as incoming instruments in the near future. 

The paper is structured as follows. In section [2] we describe 
the diff'erent models considered, in section |3] we present the re- 
sults from two diff'erent point of views: that of imaging instru- 
ments, and that of instruments with spectrometric capabilities. 
In section |4] we discuss the results and conclude in section [S] 



2. Models 

2.1. Variation of thermodynamic quantities 

In the forward modeling approach, the task of identifying MHD 
modes in the solar corona starts by implementing these modes in 
realistic numerical simulations. For the latter, it is necessary first 
to identify the theoretical observational signatures from such 
modes, free of nonlinear eff'ects. In this work we take such an 
approach and consider the simplest model of an axisymmetric 
plasma cylinder of radius R embedded in a low-jS coronal en- 
vironment, oscillating with the sausage mode, i.e. the axisym- 
metric m = mode described in [Edwin & Roberts [ ( |1983[ ). 
In this study we ignore the eff'ects of gravity, field-line curva- 
ture and twist. The variation of the thermodynamic quantities 
are found by linearizing the perturbed ideal magnetohydrody- 
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namic (MHD) equations about the magnetostatic equilibrium. In 
a cylindrical coordinate system (r, 0, z) we have: 



^ = -PoV-v 
1 

= -V/7+ -(VxB)xBo 

Ot Jd 

— = Vx(vxBo) 
ot 

dp_ ^ ypodp 

dt po dt 



(1) 
(2) 

(3) 
(4) 



where v = (Vr,V0,v^) is the perturbation velocity, p, p, B are 
the perturbed mass density, gas pressure and magnetic field, and 
Po, /?o, Bo = (0, 0, 5^,o) correspond to the unperturbed quantities 
(with [I the magnetic permeability and y = 5/3 the ratio of spe- 
cific heats). In the case of a sinusoidal perturbation leading to a 
standing mode we have V • v = -1i{r)A cosiojt) sin(fe), where A 
is the total amplitude of the perturbation and k = tiN/Lq is the 
longitudinal wavenumber (with N the longitudinal mode number 
and Lo the length of the cylinder). In the case of a sausage mode, 
"R, the radial dependence of the divergence of v must satisfy 



d^^ Id^ 2^ ^ 
ar r dr 



(5) 



Here, = 



{c]+vl){Bc^,-aj^) ' - ^s^Aiy^s + ^a)^ wnere va - 

Bz,o I ^/^P is the Alfven speed and Cs = ^|ypTp is the sound 
speed. CO satisfies the dispersion relation: 



(6) 



where ^o and Jq are the 0-th order Bessel functions. The per- 
turbed variables of interest become: 



p = sin(6i;0 sin(fe) 



P = 



AypoR . 



sm(ajt) sin(fe) 



T = sm(6i;0 sm(fe) 



co^{(jjt) sin(fe) 





AkdR 



^ cos(a;0 cos(fe). 



(7) 
(8) 
(9) 

(10) 
(11) 
(12) 



The unperturbed thermodynamic quantities are linked through 
the ideal gas law for a fully ionized hydrogen gas pq = I^pqTq. 



2.2. Base model 

Once we know the analytical solution for the thermodynamic 
quantities we produce a 3D data cube of a cylinder oscillating 
with the sausage mode. All the models considered in this work 
are variations of a 'base model'. For comparison purposes we 
choose our base model to be similar to the simulation model 
by |Gruszecki et al.] ( |2Q12| ). Accordingly, we take a discontinu- 
ous density profile with a mass density contrast Pi/pe = 5 and 



Radial veloaity 



Longitudinal velocity 




-5 
-10 
[km/s] 



[km/s] 



Density 




X 10 


-12 


1 ' 


■ 








:| 












\ \ 


1 ^ 






1 : 


1 











Temperature 




[kg r 



Fig. 1. Maps of radial velocity (top-left), longitudinal veloc- 
ity (top-right), density (bottom- left), and temperature (bottom- 
right) along t he m iddle of the cylinder for the base model - 
171 (see Sect. 2A_ for more details). The snapshot is taken at 
1/3 of the mode period. The black and white dashed lines in the 
top panels denote the 4 diff'erent line-of- sight angles considered: 
0°, 30°, 45° and 60°. For each angle two rays along the line-of- 
sight are plotted, one passing through a node of the radial veloc- 
ity (left ray of the pair) and one passing through an anti-node of 
the same quantity (right ray of the pair). The diff'erent pixel sizes 
considered are shown centred at a node: 'OR' (=25 km), 'IR' 
= 1 Mm, '2R' = 2 Mm, '3R' = 3 Mm. In the base model, IR, 
2R and 3R correspond roughly to ^/ 3, 2^/3 and i, respectively, 
where i denotes the wavelength of the mode. 



Long A mode 



Radiol velocity 



Longitudinal velocity 




Fig. 2. Same as Fig. [l]but for the long A model. See Sect, 
more details. 



2.4 



for 



a magnetic field ratio BijBe = 0.89 between the interior and 
the exterior of the unperturbed cylinder. Correspondingly, the 
Alfven speeds inside and outside the cylinder oscillate around 
350 km s"^ and 900 km s"^ respectively. The plasma p parame- 
ter is at all times lower than 0.3. The internal and external tem- 
perature for the unperturbed cylinder are set in thi s model to 
1.05 X 10^ K and 2.4 x 10^ K. For further analogy with |Gruszecki 
,et al.| ( |2012| ), we choose the wavenumber in the base model as 
kR = 2.24, where R = 1 Mm denotes the radius of the cylin- 
der. This results in a wavelength of 2.8 Mm. The amplitude A 
of the standing mode (see Eq. [7]) is chosen such as to generate 
velocity perturbations in agreement with typic al values deduced 
from spectral line broadening in the corona (Chae et al.||1998[ 
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Peter|2010| ). Accordingly, we choose A = 0.011 (similar to the 
value used in |Gruszecki et al.|2012| ). 

Typical sizes for plasma inhomogeneities in the corona are 
still not well known. This has lead to several debates concerning 
the concept of loops as monolithic or rather multi- stranded struc- 
tures in the corona. Currently, state-of-the-art coronal imaging 
devices such as SDO/AIA with instrumental point spread func- 
tions (PSF) on the order of 1 - 1.2 Mm are providing increasing 
observational evidence for a family of structures with sizes not 
resolved by the telescopes ( [Brooks et al.||2012| ). Higher resolu- 
tion imaging can be achieved with chromospheric instruments 
through tracing of coronal rain in loops, which has provided evi- 
dence for strand like structures on the order of a few hundred km 
( Antolin & Rouppe van der Voort|2012 (). While these works re- 
veal sub-arcsecond density structure providing support for multi- 
stranded loop models, it is reasonable to think that the collection 
of strands composing the loop (defining the main waveguide) 
is actually much larger and behave collectively, agreeing with 
our choice of 1 Mm for the radius of our loop. Such families of 
neighbour strands traced by coronal rain indeed h ave s ev eral Mm 
in diameters, as seen in Antolin & Verwichte (2011) o jAntolin] 
& Rouppe van der Vo ort (2012 ). Furthermore, it has been shown 
that the characteristics of a sausage mode are rather insen sitive to 
the presence of cross- sectional inhomogeneities in loops ( |Pascoe| 
let al.|2007al [2009 ^. 

Numerically, our numerical model consists of a 204x204xA^^ 
grid. A^^ denotes the number of points in the z direction and 
is chosen so that the diff'erent line-of-sight rays for the inten- 
sity calculation (explained in Sect |2.3| ) traverse the same cross- 
sectional area in the model. In the base model, three wavelengths 
are enough, leading to = 342 points. The spatial resolution 
in X, y and z axes is the same: about 25 km. The fast sausage 
mode in the corona produces very low cross-sectional area vari- 
ation, which for some of the considered models leads to very 
low intensity variations. To correctly calculate the variations in 
those models we have also considered a numerical model with 
10 times the resolution in each direction. Correspondingly, all 
results related to intensity modulation are obtained with the high 
resolution models. All other results do not vary significantly be- 
tween the low and high spatial resolution models. 

Figure [T] shows a snapshot of our base model, representing 
a cut along the middle axis of the cylinder. The sausage mode 
creates in-phase perturbations between the density and the mag- 
netic field, and anti-phase perturbations between the density and 
the radius. Furthermore, the radial and longitudinal velocities are 
71 / 4 out of phase with the density. The snapshot in Fig.[T]is taken 
at a time of 1/3 of the period P, which, in this model, corresponds 
to t ^ 4.6 s. Positive velocity values correspond to outflow to- 
wards the axis of the cylinder. From Fig. [T] it can be noticed 
that the perturbation of longitudinal velocity is about 1 order of 
magnitude smaller than the radial perturbation. They have max- 
imum values of about 19 km s~^ and 3.4 km s~\ respectively. 
At the location of an anti-node of the radial velocity, density os- 
cillates about 7 % around a mean value of 4.9 x 10"^^ kg m"^, 
and temperature oscillates about 5 % around a mean value of 
1.048 xlO^ K. 



2.3. Calculation of synthetic spectra 

For the calculation of synthetic spectra we choose difl'erent view- 
ing angles on the structure, simulating difl'erent line-of-sight 
configurations. Denoting by 6 the angle between the line-of-sight 
and the perpendicular to the axis of the cylinder, we consider 
the cases: ^ = 0°, 30°, 45°, 60°, shown in white or black dashed 



lines in Fig. [T] For each line-of-sight, we define ray 1 and ray 
2 as two particular rays crossing the cylinder at a node and an 
anti-node of the radial velocity, respectively (unless stated oth- 
erwise, the terms 'node' and 'anti-node' will refer to the radial 
velocity amplitude throughout the manuscript). For simulating 
observations it is futher important to consider the spatial, tem- 
poral and spectral resolution of the instrument at hand. For each 
ray we consider 4 diff'erent resolution 'pixels': 'OR' (0 radius, 
equal to the spatial resolution of our model, i.e. 25 km), 'IR' (1 
radius ^ 1/3^, where £ is the wavelength of the mode in the base 
model), '2R' (2 radii ^ 2/3^) and '3R' (^ 0- These pixel sizes 
correspond to the dotted squares in Fig.[T] A pixel of 1 R would 
correspond roughly to the size of the instrumental PSF for the 
171 and 193 channels of the AIA instrument of SDO ( |Grigris| 
et al.|2012, AIA Team Online Document| ). Taking into account 
that the timescale of the sausage mode in a coronal structure such 
as this one is on the order of seconds (the period P of the mode 
in the base model is about 5 s), which is very low with respect 
to the current cadence of available instruments (at most 1 s for 
imaging instruments, in general 12 s), we consider the follow- 
ing temporal resolutions: 'Os' (equal to the temporal resolution 
in our model, i.e. 0.15 s), '1 s' (^ P/5) and '2 s' (^ 3P/7). The 
spectral resolution in our numerical model is set at 1 .4 mA for 
the Feix 171 line and 1 mA for the Fexii 193 line, leading to a 
spectral sampling of 48 points across 99 % of the spectrum (de- 
tails about the selected lines are given below). We also consider 
coarser spectral resolutions of A A = 4.2 mA and A A = 8.5 mA, 
leading to samplings of 16 and 8 points respectively. 

The quantity that we are mostly interested in synthesising is 
the specific intensity for optically thin (coronal) plasma along a 
given line-of-sight: 



(13) 



where Ab is the (coronal) abundance of the element (denoted by 
its wavelength at line centre Aq) relative to hydrogen, geff is the 
Gaunt factor, is the oscillator strength of the line, G^^^iT, Ue) 
is a mostly temperature dependent function characterising the 
response of the plasma for that particular line, tie is the elec- 
tron number density, and where the integral is evaluated along 
a ray crossing the emitting plasma {dl is an element of length 
along the line-of-sight). The G^^iT, Ue) function depends mainly 
on the ion population fraction, which can be obtained by solving 
the time-dependent ion balance equations. While such a treat- 
ment would provide a thorough analysis of the radiative response 
of the plasma, it would take us outside the scope of this paper 
(such a treatment has been followed by, e.g. Bradshaw & Mason 
2003 ). Here, we will assume that the plasma is at all times in 
equilibrium ionisation and we will estimate the eff'ects of non- 
equilibrium on the observable quantities. 

The purpose of this work is not to determine exactly the ab- 
solute value of line intensities, but rather to analyze the eff'ect on 
the intensities from temperature, density and especially the ge- 
ometry (the line-of-sight rays). Hence, we concentrate here on 
the quantities that compose the integrand in Eq. 13 The calcu- 
lation of the integral in Eq. [13] for all rays crossing the domain 
is a computationally expensive procedure that can be consider- 
ably sped up with the following techniques. For the calculation 
of Gxq(T, He) we produce a look-up table depending o n tempera- 



ture and density using the CHIANTI atomic database ( Dere et al. 
1997 , 2009). For the latter we choose the chianti collisional ioni- 
sation equilibria file (chianti . i oneq) and coron al abundances 
(sun_coronal . abund file, from Feldman 1992| ). The look-up 
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Base model - 171 

l.o.s. integrated intensity 0° case l.o.s. integrated intensity 30° case 





2 4 6 

(perp. to l.o.s)/R 



l.o.s. integrated intensity 45° cose l.o.s. integrated intensity 60° case 
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Base nnodel - 193 

l.o.s. integrated intensity 0° case l.o.s. integrated intensity 30° case 





l.o.s. integrated intensity 45° cose l.o.s. integrated intensity 60° cose 




4 6 
(perp. to l.o.s)/R 




Fig. 3. Integrated intensity along the line-of-sight for the base 
model - 171, at the 4 different angles considered: 0°,30°,45° 
and 60° (see Sect. |2.4| for more details). The snapshot is taken 
at 1/3 of the mode period. The crosses in each pan el de note the 
location of the 2 line-of-sight rays defined in Sect. 2.3 In each 
panel, the left and right crosses correspond to the rays crossing 
a node and an anti-node of the oscillation. 



table is composed by a 2000 x 100 points uniform sampling in 
the electron number density Ue - temperature T domain (respec- 
tively), where logn^ g [8,10] and logT g [4.95,6.5], which 
covers the values obtained in our numerical model. Since the 
G function is mostly dependent upon temperature the values of 
the table are only interpolated to the temperature values of each 
pixel in our model. We produce histograms of the line-of-sight 
velocity choosing a bin size of 0.5 km s"^ (much smaller than the 
plasma motions produced by the perturbation) and work with the 
resulting bins instead of pixel-by -pixel. For each resulting bin 
we produce a second histogram of the emissivity G^^iJ ,ne)n^^ 
with the data points within the bin. Thanks to the symmetry in 
our model the bin number necessary for high accuracy for this 
second histogram is lower than 10000. The spectral line is then 
calculated for each emissivity bin by assuming a Gaussian dis- 
tribution for the plasma. The maximum of the Gaussian is given 
by the value of the emissivity bin. The width of the Gaussian 
is set by the thermal broadening, given by the logarithmic av- 
erage temperature of the pixels within the emissivity bin. The 
displacement of the Gaussian's centre is set by the bin value of 
the velocity histogram. The last step consists in integrating along 
the line-of-sight the Doppler shifted emission. This involves an 
interpolation from the values of the cartesian grid of the Doppler 
shifted emission to equally spaced points along the line-of-sight 
rays. With the outlined method the calculation time of the syn- 
thetic spectrum for one line for 4 different line-of-sight angles is 
8 to 10 times faster than with a normal pixel-by-pixel method. 

2 A. Models 

The modulation of the intensity due to the sausage mode is anal- 
ysed in several models which differ in 1 specific way from our 
base model. 

• The 'base model - IJV - We analyze the modulation in in- 
tensity for our base model in the transition region line Fe ix 
171.073 A, which has a formation temperature at logT = 
5.9 K. 



2.4 



Fig. 4. Same as Fig. ijbut for the base model - 193 (see Sect, 
for more details). The pixel sizes corresponding to IR, 2R and 
3R are shown around the location of ray 2 (anti-node). 

• The 'base model - 193^ - We analyze the modulation 
in intensity for our base model in the coronal line Fexii 
193.509 A, which has a formation temperature at logT = 
6.1 K. 

• The 'high-T modeV - The effects of high external temper- 
ature are analysed in this model. We use an external tem- 
perature of 8 X 10^ K, which results in a slightly higher in- 
ternal density, with a mean value at the anti-node of 5.5 x 
10"^^ kg m" and a variation around this value of roughly 
7%. 

• The 'long A modeV - Keeping the same values of internal 
and external thermodynamic quantities for the unperturbed 
cylinder (base model), here we model the longest possi- 
ble wavelength for the sausage mode for these conditions: 

5 Mm, which results in a wavenumber of kR = 1.25. 
Figure |2] shows a snapshot of a cut along the middle axis of 
the cylinder for the radial and longitudinal velocities, den- 
sity and temperature for this case. The snapshot is taken at 
1/3 of the period, which in this case is 5.8 s. The white and 
black lines show the line-of-sight rays at the different angles 
considered in this study. For this model two wavelengths in 
the z direction are enough to cover the line-of-sight rays, giv- 
ing A^^ = 408 grid points. The density and temperature have 
roughly the same average values as before but now exhibit 
slightly larger variations at the anti-nodes, about 8.5% and 
6% respectively. The radial and longitudinal velocities have 
now maxima at 23 km s~^ and 3 km s~^ . 



3. Results 

Since the aim of this work is to identify the observational signa- 
tures of the sausage mode in coronal structures we present the 
results of our models as they would be perceived by the instru- 
ments at hand. We therefore divide the results into 2 relevant 
categories depending on the instrument: imaging and spectro- 
scopic. 

3.1. For imaging instruments 

Figure [3] shows the integrated intensity for the base model - 171 
(corresponding to the 10 times higher spatial resolution model), 
obtained by integrating the intensity in Eq. [13] over the full Fe ix 



5 



Antolin & Van Doorsselaere: Intensity perturbations by sausage modes 




Fig. 5. Integrated intensity along the line-of-sight (normalised with respect to the equilibrium value) for the 0° angle case for the 
base model - 171 (left panels) and the base model - 193 (right panels) for two specific rays: ray 1 crossing a node (top most panels) 
and ray 2 crossing an anti-node (bottom panels), shown with crosses in Figs. [3] and |4] Diff'erent colors denote diff'erent pixel sizes: 
OR, IR, 2R and 3R in red, purple, blue and green, respectively, where R = 1 Mm ^ £ denotes the radius of the cylinder and OR 
denotes the spatial resolution of our model: 25 km. The black curve denotes the density profile at the centre of the cylinder and the 
same location along the axis of the cylinder as rays 1 and 2. The intensity range of the top panels is set significantly smaller than 
that of the bottom panels in order to show the low intensity modulation at the nodes. 



171 A line, for the 4 line-of-sight angles: 0°,30°,45° and 60°. 
The crosses in the centre of each panel show the location of rays 
1 and 2 crossing the cylinder defined in Sect. |2.3[ The snapshot 
is taken at 1/3 of the period, coinciding with that of Fig. [T] 

One of the main features of this figure is the absence of sig- 
nificant intensity structuring along the cylinder, in sharp contrast 
with the density or temperature maps shown in Fig.[T] The inten- 
sity pattern matches the radial expansion of the cylinder, which 
is barely discernible, equal to only 1 pixel in the low resolution 
model (25 km). We thus have a longitudinal structuring in in- 
tensity qualitatively similar to the longitudinal structure of the 
radial velocity. On the other hand, contrary to what would be 
expected if the intensity is assumed to be proportional to the 
density squared, the regions of increased intensity correspond to 
regions of low density at the centre of the cylinder. This is due 
to the fact that the Feix 171 A line has a formation temperature 
at logT = 5.9, and the sausage mode is only slightly compress- 
ible in the corona, producing temperatures in our cylinder that 
oscillate between logT = 5.99 and logT = 6.04. Therefore, 
a drop in temperature in our model leads to an increase in the 
fraction of the Feix population, thus increasing the emission, 
as can be seen in Figure [6] This eff'ect is further reinforced by 
the fact that at 0° viewing angle the rays crossing these regions 
will cross more material since they correspond only to regions 
of expansion (lower temperature). On the other hand, as shown 
in Figure [6j the Fexii 193 A line has a formation temperature 
at logT = 6.2 K. In this case, higher temperature implies an 
increase in the ion population, leading to higher intensity and 
thus creating a structure that is similar to the temperature or den- 
sity structure, as can be seen in the base model - 193 shown by 
Figure |4] This fact highlights the dominant role of the contri- 



bution function G^^iT^ne) in the formation of the emergent in- 
tensity. It is important to remind the reader that these results are 
subject to the assumption of equilibrium ionisation, and there- 
fore that the plasma responds instantaneously to the changes in 
temperature and density from the MHD mode. The limitations 
of this assumption will be discussed at the end of this section. 

Figure[5]shows the evolution of the integrated intensity along 
the line-of-sight at an angle of 0° for the base model - 1 71 (left 
panels) and the base model - 193 (right panels), for the rays 1 
(top panels) and 2 (bottom panels) crossing at a node and an anti- 
node of the oscillation, respectively. Any other location along 
the axis of the cylinder displays features similar to either ray 1 
or ray 2. For each ray diff'erent spatial resolution pixels are taken, 
marked in diff'erent colors (see figure for details). In the case of 
the base model - 171 have a very low modulation of the in- 
tensity, as compared to the density modulation at the centre of 
the cylinder (black curve in the figure). While the latter varies 
by < 7%, the intensity varies by < 2% at an anti-node and is ba- 
sically negligible at the location of a node. Furthermore, as ex- 
plained earlier, the variation of the intensity is in antiphase with 
that of the density due to the fact that the contribution function 
Gi7i(r, He) has a peak at lower temperature values than those of 
the plasma in the cylinder. The slightly longer column thickness 
created by the area expansion creates the slight asymmetry be- 
tween the two dips of the intensity cycle. 

This scenario is in sharp contrast with that of the base model 
- 193, shown in the right panels of Fig. [5] The Gig^iT, n^) func- 
tion presents a peak at a higher temperature value than that of the 
plasma inside the cylinder. In this case we thus have in-phase in- 
tensity variations with the density (or temperature) profile. Also, 
as shown in Figure[6j the GigsiT, rie) contribution function is sig- 
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G,7i(T.ne) 




Fig. 6. Contribution function G^Q(r,^g) for Feix 171 and 
Fexii 193. The maximum formation temperature for these ions 
are, respectively, log T = 5.9 and log T = 6.2 K. The black paral- 
lelogram in each surface has its corners defined by the minimum 
and maximum temperature and density values found in the base 
model, and therefore includes the G^o(r,^e) values used here. 
Note that in order to show the parallelogram for the Gi7i(r, rie) 
surface the temperature range has been inverted. 



nificantly more sensitive to the temperature changes inside the 
cylinder than the Gi7i(r, rie) function, leading to intensity varia- 
tions of ^ 35 %. The diff'erence in intensity from the initial value 
is greater for the maximum values than the minimum values. 
This is due to the fact that the gradient of the G^^Cr, rie) function 
is steeper for both 171 and 193 above its initial value (mainly in 
the direction of higher temperature) than below the initial value. 
This eff'ect is further amplified by the density enhancement at a 
time of contraction, since each pixel along the line-of-sight con- 
tributes with the squared of the local density value. We further 
notice 2 small bumps for ray 1 at times of density extrema when 
considering pixel sizes diff'erent from OR. These are again due 
to the combined eff'ect of the steepness of the G^^Cr, rie) gradient 
and enhanced density in regions of contraction, which are within 
the line-of-sight for coarse pixel sizes. 

The variation of the intensity decreases significantly with the 
pixel size for all models, especially at anti-node locations when 
viewing at a 0° line-of-sight angle. In the considered models the 
intensity variation is reduced by factors of 10 - 17 when increas- 
ing the pixel size to 3R. This is expected since a larger pixel 
size implies a larger probability of crossing mixed regions of 
expansion and contraction. For the models with a wavenumber 
ka = 2.24, the unit of the pixel is slightly larger than | of the 
wavelength, and therefore 3 pixel units 3R) not only covers 
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Fig. 7. Comparison between an intensity calculated as / p^l, 
where / is the column thickness along the line-of-sight (solid), 
the base model - 171 intensity hnit) (dotted) and the base model 
- 193 intensity hgsit) (dashed) for a viewing angle of 0° and a 
pixel size resolution of OR for ray 2. The diff'erences between the 
lines illustrates the importance of the Gx^iT, rie) function. 



Integrated intensity along l.o.s. OR pix 



Base model - 1 93 




Time [s] 

Fig. 8. Same as the right panel of Fig.jijbut fixing the pixel size 
resolution to OR and varying the angle of the line-of-sight for 
the base model - 193. 0°, 30°, 45° and 60° correspond to the red, 
purple, blue and green curves, respectively. 



the entire region of expansion and contraction, but also produces 
a slight imbalance between the two, as can be seen in Figs.|3]and 
[4] Hence, at a time of contraction the pixel will cover slightly 
more expanding material, and vice- versa, resulting in a reversal 
of the variation in intensity at anti-node locations (ray 2). This 
can be seen in Fig.|5]for both models, 171 and 193. However, in 
this case, most of the longitudinal intensity structuring is erased, 
and with it the diff'erence between the 171 and the 193 line in- 
tensity variation. Furthermore, if our pixel is centred in a node, 
regardless of the size of the pixel the amount of expansion and 
contraction will always be similar, thus producing minimal in- 
tensity variation. As can be seen in the top right panel of Fig. [5] 
this variation attains a maximum for a 2R pixel size, which is 
linked to the fact that the amount of material entering and leav- 
ing the line-of-sight is close to a maximum for this size. This is 
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of course dependent on the shape of the pixel, which in our case 
is a square. 



0° case 

Velocity [km/s] Velocity [km/s] 

-40 -20 20 40 -40 -20 20 40 




171.06 171.08 171.1 171.06 171.08 171.1 

Wavelength [A] Wavelength [A] 

Base nnodel - 171 



0° case 

Velocity [km/s] 2R res Velocity [km/s] 
-40 -20 20 40 -40 -20 20 40 




Fig. 9. Time variation of the specific intensity h integrated along 
ray 1 (left panels) and ray 2 (right panels) for the base model - 
1 71 when viewing at an angle of 0° and for a pixel size of OR 
(upper panels) and 2R (lower panels). Contours of the intensity 
are also shown for better visualization of the variation. 



As explained in the introduction, a common approximation 
to intensity calculation is / p^/, where / is the column thickness 
along the line-of- sight. In Fig. [7] we compare this approach to 
the outcomes of the more proper intensity calculation performed 
in this work, plotting p^l and the results from the base model 
- 171 and base model - 193. To clearly show the diff'erence we 
plot the case of a 0° viewing angle and a pixel size resolution 
of OR for ray 2 (anti-node). By following the I ^ p^l approxi- 
mation we find maximum intensity variations below 3 %, which 
match those obtained with the 171 line. However, both profiles 
are in antiphase with each other due to the fact that the Fe ix ion 
population decreases with increasing temperature for the range 
of temperatures considered here. The discrepancy in amplitude 
with the 193 line intensity is obvious. This illustrates the impor- 
tant role of the contribution function G^Q(r, rie) in the formation 
of the emergent intensity. 



45° case 

Velocity [km/s] 0^ Velocity [km/s] 

-40 -20 20 40 -40 -20 20 40 




45° case 

Velocity [km/s] 2R res Velocity [km/s] 
-40 -20 20 40 -40 -20 20 40 



4 



3 

.§ 2 




171.06 171.08 171.1 171.06 171.08 171.1 

Wavelength [A] Wavelength [A] 

Base model - 171 



Fig. 10. Same as Fig.[9|but for a line-of-sight crossing the cylin- 
der at an angle of 45°. 



The variation in the integrated intensity due to diff'erent line- 
of-sight angles is shown in Fig. [8] In order to show best the vari- 
ation we fix the pixel size resolution at OR and consider the base 
model - 193, for which large intensity variation is found for the 
temperature and density variations created by the sausage mode. 
By increasing the line-of-sight angle the rays cross more regions 
of hot (contraction) and cold (expansion) plasma, thus reducing 
the variation of the intensity by a factor of 2.5 (down to ~ 10 %). 
For rays crossing nodes (ray 1), as discussed previously, the vari- 
ation is greatly reduced from the crossing of mixed temperature 
regions in equal quantities along the rays. As in Fig. [5] (right 
panel), two small bumps can be seen for non-zero viewing angles 
produced by the combined eff'ect of the local gradient steepness 
of the G^Q(r, ^g) function and the regions of enhanced density 
along the line-of- sights. 

Apart from the 2 models discussed in this section, the other 
models exhibit the same behaviour as that of the base model - 
171, with minor quantitative diff'erences. The main diff'erences 
are discussed in section [34l 
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Fig. 11. Left panel: Pixel size resolution dependence of the time variation of the specific intensity h for the 0° angle case of the 
base model -171. From top to bottom we show quantities ensuing from a Gaussian fit to the spectral line. The first 3 upper panels 
correspond to the Doppler velocity, the line width and the goodness-of-fit parameter for the case of ray 1 (node). The 3 lower 
panels show the same quantities for the case of ray 2 (anti-node). Diff'erent line styles correspond to diff'erent pixel size resolution: 
OR (solid), IR (dotted), 2R (dashed), 3R (dot-dashed). The line width is normalised with respect to the initial value at time ^ = 0. 
The unit for thex^ parameter is the logarithm of the normalised value, where the normalisation factor is the equilibrium value. Right 
panel: Angle dependence of the time variation of the specific intensity for the OR pixel size case of the base model -171. The 
arrangement of these panels is the same as for the left panels. Diff'erent line styles denote diff'erent line-of-sight angles: 0° (solid), 
30° (dotted), 45° (dashed), 60° (dot-dashed). For both left and right panels the variation in time of density (red curve) and radial 
velocity (blue curve) at the centre of the tube is plotted for reference in arbitrary units in the Doppler velocity panel. 



Non-equilibrium ionisation effects 

Perhaps the most important assumption in the intensity modula- 
tion presented here is to consider that the plasma is at all times 
in equilibrium ionisation, i.e. that the plasma response to the 
temperature and density changes produced by the MHD mode 
is instantaneous. However, this is never the case and significant 
departures from equilibrium can occur in timescales which can 
be short with respect to the timescale of the driver (|De Moortell 
& Bradshaw||2008| ). The population densities of the ions in the 



plasma will vary according to the ionisation and recombination 
rates, which depend on the thermodynamic conditions at each 
specific moment of time. The rigorous solution to this prob- 
lem involves solving the time-dependent ion balance equations, 
which would take us outside the scope of this paper. However, 
we can estimate the consequences of non-equilibrium ionisation 
on the emergent intensity. For this we assume that the plasma 
follows a Maxwellian distribution at all times. Timescales for 
ionisation and recombination equilibrium for the ions Fe ix and 
FeXII can then be obtained through CHIANTI. These values 
are shown i n Table p] (in the case of a /c-distribution, as shown by 
[Dzifcakova 2002[ , the rates would increase leading to smaller 
timescales than those presented here). 

For Fe ix we can see that the ionisation and recombination 
timescales to reach equilibrium are on the order of 2 s, which is 
almost half the period of the sausage mode. This means that at 
a time of expansion (cooling) the material will just have started 
being ionised and will therefore start recombining without at- 
taining a full ionisation. Hence, we have a time lag eff'ect and an 



attenuation (and deformation) of the intensity profile. The loss in 
amplitude will be similar to that obtained when integrating the 
profile over a period of 2 s, in which case the intensity is reduced 
by a factor of roughly 1.5. This eff'ect is shown in Figure 13 For 
Fe XII the ionisation and recombination timescales are between 
13 s and 20 s, that is, considerably larger than the period of our 
mode. In this case, there would be basically no modulation of 
intensity. 

Non-equilibrium ionisation eff'ects would however have no 
influence on the spectrometric results presented in section |3]2 
Indeed, as we will see in the next section, most of the spectra 
broadening is of turbulent nature and not of thermal origin. 

3.2. For spectroscopic instruments 

We now analyze the spectral variation caused by the sausage 
mode, which is of interest for observations performed by spec- 
troscopic instruments. For this we concentrate on rays 1 and 2 
only, simulating slits that would cross the cylinder at a perpen- 
dicular angle, centred at a node or an anti-node, respectively. All 
other positions along the cylinder present features similar to ei- 
ther ray 1 or ray 2. 

Figures [9] and 10 show the time variation of the specific in- 
tensity for the base model - 171 for the line-of-sight angles of 
0° and 45° respectively. Upper and lower panels in the figures de- 
note, respectively, OR and 2R pixel size resolution, and left and 
right panels denote, respectively, ray 1 (node) and ray 2 (anti- 
node). For the 0° angle and 2R pixel size case we have periodic 
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Table 1. Estimates of timescales for ionisation (I) and recombination (R) equilibrium for Fe ix and Fe xii inside the model 



Ion 


T^,n ^ 1 X 10^ K 

He, min ^ 2.7 X 10^ Cm"^ 


T^ax ^ 1.1 X 10^ K 
ne, ma.. ~ 3 X 10^ cm"^ 


Fe IX 
Fe XII 


2s(I), 2.1 s(R) 
22.5 s (I), 1.5 s (R) 


1.3 s (I), 2.1 s(R) 
12.9 s (I), 1.5 s (R) 



Notes. The timescales are given in seconds and are calculated as l/(We x X{T)), where X(T) denotes either the ionisation or recombination rate 
obtained with CHIANTI assuming a Maxwellian distribution for the plasma. The temperature (T) and number density (Ue) values used for the 
calculation are extremum values found in the centre of the cylinder. 
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Fig. 12. Estimation of the spectral line broadening contribution 
from the unresolved velocity perturbations (which we denote as 
turbulence) and temperature perturbations. The Gaussian width 
corresponding to the thermal broadening contribution (dashed 
curve) is obtained from the average of temperature for the pix- 
els along the line-of-sight crossing the cylinder. The Gaussian 
width corresponding to the turbulent broadening contribution 
(solid lines) is obtained from a Gaussian fit to the sum of multi- 
ple Gaussians, one for each pixel in the cylinder along the line- 
of-sight, with a fixed temporal averaged thermal broadening and 
shifted centre given by the line-of-sight velocity at that pixel. 



spectral broadening for rays crossing the cylinder at any loca- 
tion, especially clear for the OR pixel size case of ray 2. This fea- 
ture is only absent for the case of ray 1 with OR pixel size. This 
is expected since for the latter the pixel size does not include 
any ray crossing a region of expansion or contraction, contrary 
to all the other cases. When watching at an angle, rays will cross 
regions of expansion and contraction by diff'erent amounts when 
the slit is centred at a node, and in equal amounts when centred 
at an anti-node. This creates an asymmetry for the former and 
therefore an overall dominating redshift or blueshift. This trans- 
lates in clear periodic blueshift or redshift excursions for the OR 
pixel size case of ray 1 (left upper panel of Fig. 10), but hardly 
discernible for the 2R pixel size case (lower left panel). As ex- 
pected, the periodic broadening is also present for this viewing 
angle for any pixel size. 

The pixel size resolution dependence of the time variation 
of the specific intensity observed in Figs. [9] and [T0| is shown 
in the left panels of Fig. [TT] for the 0° angle case. In these pan- 



els we show the result of Gaussian fits to the specific intensity 
7^(0 for all the considered pixel sizes. The time variation of the 
Doppler velocity, width of the Gaussian fit, and the goodness- 
of-fit statistic are shown from top to bottom in the first 3 
panels for ray 1 and the last 3 panels for ray 2. Since the 0° 
angle case is considered and the sausage mode is mainly a radial 
mode, negligible Doppler velocity results from the Gaussian fit 
for any pixel size resolution. On the other hand, the line width 
and the non-Gaussianity parameter present significant varia- 
tion, and have profiles that are in phase with the square of the 
radial velocity. This means that they both have a double peri- 
odicity. Maximum variations close to 20 % and 50 % are ob- 
tained for the line width, respectively for rays 1 and 2. Notice 
that for the case of ray 1 (ray 2) the spectral line broadening 
increases (decreases) for pixel sizes up to 2R. For 3R, which 
corresponds to 1 wavelength, and therefore an equal amount of 
expanding and contracting gas, the broadening is reduced (in- 
creased) for ray 1 (ray 2). The decreasing factors range between 
2.5 and 3.4 between the models. The temperature (or density) 
modulation by the sausage mode contributes to a certain extent 
to the non-Gaussianity of the line profile. This is especially the 
case for large pixel sizes (> 2R), for which the contribution from 
the radial velocity is lowest due to the large range of velocities 
included, producing a broad Gaussian-like profile. This can be 
seen in the;^^ profiles as a bump at times of maximum and min- 
imum densities at the centre of the cylinder. These bumps are of 
the same order as the perturbations created by the radial veloc- 
ities, and therefore can give to an observer the impression of a 
quadruple periodicity for the;^^ parameter (period of P/4 where 
P is the period of the sausage mode). 

The angle dependence of the time variation of the specific 
intensity for rays 1 and 2 is shown in the right panels of 
Fig. [TT] for the OR pixel size resolution case of the base model 
-171. For non-zero viewing angles the rays now cross both re- 
gions of expansion and contraction. The Doppler velocity values 
ensuing from the Gaussian fits are therefore the emission mea- 
sure weighted averages over these regions. A small contribution 
comes also from the longitudinal velocity component, but it is 
however one order of magnitude smaller than the transverse ve- 
locity component. As observed in Fig. 10 depending on the line- 
of-sight angle one can have at a given time either a dominating 
blueshift (30° and 45°) or redshift (60°). The variations in spec- 
tral line broadening given by the line width and the;^^ parameter 
decrease significantly with the angle when crossing an anti-node 
(with maximum factors for line width between 5.6 and 11.3 be- 
tween the models). On the other hand, rays crossing a node have 
slightly increasing variations of spectral line broadening with an- 
gle (below 10 %). The maximum velocity component is directed 
perpendicularly to the cylinder axis. The components creating 
the width of the spectral line therefore decrease with the sine of 
the line-of-sight angle. The contribution to the non-Gaussianity 
of the line from the temperature and density profiles can also be 
observed at various angles. The bumps in the;^^ parameter pro- 
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file created at times of maximum or minimum temperature (and 
density) at the centre of the cylinder are of the same order of the 
perturbation created by the radial motions for any angle in the 
case of ray 1 and for angles > 45° for ray 2. 

In order to estimate the respective contribution from turbu- 
lence (unresolved velocity perturbations) and temperature to the 
broadening, we plot in Figure[T2]the time variation of 2 Gaussian 
widths obtained from rays 1 and 2 with a viewing angle of 0° and 
and pixel size of OR. The Gaussian width corresponding to the 
thermal broadening contribution (dashed curve) is obtained from 
the average of temperature for the pixels along the line-of-sight 
crossing the cylinder. The Gaussian width corresponding to the 
turbulent broadening contribution (solid lines) is obtained from 
a Gaussian fit to the sum of multiple Gaussians, one for each 
pixel in the cylinder along the line-of-sight, with a fixed tempo- 
ral averaged thermal broadening and shifted centre given by the 
line-of-sight velocity at that pixel. We can see that the turbulent 
broadening is the main contributor to the spectral line broaden- 
ing. Only for cases of little plasma motion along the line-of-sight 
(0° angle and OR pixel size, or at times of radial velocity) is the 
thermal broadening comparable to the turbulent broadening. 
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3.3. Effects of temporal and wavelength resolutions 

In section [TT] we discussed the eff'ects from non-ionisation equi- 
librium, and in the last section we showed the eff'ects of spatial 
resolution and line-of-sight angle on the detection of intensity 
modulation by the sausage mode. In this section we would like 
to present the main eff'ects due to temporal and wavelength res- 
olutions. The temporal and wavelength resolutions of our base 
model are, respectively, 0.15 s and 1.4 mA. Here we will discuss 
the diff'erences introduced by taking temporal resolutions of 1 s 
and 2 s (0.15 s, 1 s and 2 s correspond, respectively, to 3 %, 
22 % and 44 % of the mode period), and wavelength resolutions 
of 4.2 mA and 8.5 mA (1.4 mA, 4.2 mA and 8.5 mA correspond, 
respectively, to 48, 16 and 8 points sampling of the spectral line). 



Figure 13 shows the eff'ect of temporal resolution on the inte- 
grated intensity and spectral broadening obtained from the mod- 
ulation by the sausage mode for the base model -193. We choose 
a viewing angle of 0° and ray 2 (anti-node), for which the ef- 
fects are strongest for the integrated intensity (top panel), line 
width (third panel from top to bottom), and parameter (bot- 
tom panel), and a viewing angle of 45° and ray 1 (node) for the 
Doppler velocity (second panel from top to bottom). We also 
fix the spatial resolution to OR. Diff'erent temporal resolutions 
are obtained by integrating the intensity profiles over temporal 
windows corresponding to the designated resolution. Since the 
period of the sausage mode is about 4.5 s only a few points of 
the full profile are obtained per period. In order to obtain full 
period profiles approximately 7 and 13 periods are needed, re- 
spectively, for 1 s and 2 s resolutions. Longer time integrations 
also generate phase diff'erences between the profiles equal to the 
resolutions of each profile. We choose however to correct for 
this eff'ect and present the shifted profiles in Fig. [13] in order to 
facilitate comparison. 

the intensity varia- 



As seen in the top panel of Fig. 13 



tion from the sausage mode decreases with temporal resolution. 
Factors ranging from 1.2 to 1.5 were found among the models. 
Although this eff'ect may not appear important for strong inten- 
sity variations as those obtained for the base model - 193, they 
are especially relevant for low intensity variations, as those ob- 
tained in the base model -171. Longitudinal variation can there- 
fore be erased in this case. 



Fig. 13. Eff'ects of the temporal resolution on the modulation of 
observable quantities. We take the base model - 193, fix the pixel 
size resolution to OR and consider cadences of 'Os' (= 0.15 s, 
solid curve), 'Is' (= 1 s) ^ P/5 (dotted curve) and '2s' (= 2 s) ^ 
3P/7 (dashed curve), where P ^ 4.5 s is the period of the mode. 
The upper panel shows the temporal evolution of the normalised 
integrated intensity along the line-of-sight for ray 2 taking a 
viewing angle of 0°. The rest of the panels, respectively from top 
to bottom, show the temporal evolution of the Doppler velocity, 
the line width and (logarithm of) the;^^ non-Gaussianity param- 
eter resulting from Gaussian fits to the specific intensity profile 
7^(0 for the diff'erent temporal resolutions. In order to show best 
the eff'ects on the Doppler velocity, only for this quantity we take 
the case of ray 1 with a viewing angle of 45°. The normalisation 
factor for the line width is the initial value at time t=0 for the Os 
temporal resolution. That of the;^^ parameter is the time average 
value of the profile for ray 1 for Os temporal resolution. 



In the lower 3 panels of Fig. 13 we show the eff'ects of tem- 



poral resolution on the Doppler velocity, the line width and the 
non-Gaussianity parameter obtained from the Gaussian fits 
to the specific intensity profiles hit). The Doppler velocity, be- 
ing the displacement of the bulk of the plasma along the line- 
of-sight, is the least aff'ected quantity, with maxima decreasing 
from 7 km s"^ to < 5 km s"^ (decreasing factors between 1.1 to 
1.5 among the models). On the other hand, the strength of the 
spectral broadening depends mostly on the fast components of 
the plasma. Correspondingly, the line width variations of about 
40% obtained in the model are reduced to 30% and < 10%, re- 
spectively, for time cadences of 1 s and 2 s, compromising its 
detectability with the current instruments. Accordingly, decreas- 
ing factors between 4.6 and 6.3 are found between the models. 
Similarly, the variation of the x^ parameter is reduced (spectra 
become more and less Gaussian at times of expansion and com- 
pression, respectively), especially in the 2 s cadence case, where 
it is reduced by several orders of magnitude. Hence, for the 2 s 
time cadence, the spectral line is seen as a profile with almost 
constant shape and width roughly equal to the average width of 
the true profile, but for which bulk Doppler motions can still be 
observed. 
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Fig. 15. Model comparison of the temporal variation of the spectral line broadening. The line width (upper panels) and (logarithm 
of) the non-Gaussianity parameter (lower panels) resulting from the Gaussian fits to the specific intensity 7^(0 for the OR pixel 
size case are presented. From left to right the first two panels correspond to rays 1 and 2 for the 0° angle case. The last two panels 
correspond to rays 1 and 2 for the 45° angle case. Diff'erent line styles denote the diff'erent models studied in this work, as indicated 
in the legend. The normalisation values are the time averages of the corresponding quantities for ray 1 and 0° viewing angle. 



:ose model - 193 Wavelength dependence 




Fig. 14. Temporal evolution of the Doppler velocity (upper 
panel), the line width (centre panel) and (logarithm of) the 
non-Gaussianity parameter (lower panel) resulting from 
Gaussian fits to the specific intensity profile lx(t) obtained from 
the base model -193, for diff'erent wavelength resolutions. We 
consider samplings of 1.4 mA (solid curve), 4.2 mA (dotted 
curve) and 8.5 mA, leading to 48, 16 and 8 points along the 
spectral line, respectively. In order to show best the eff'ects of 
the wavelength we fix the pixel size resolution to OR and show 
the case of ray 1 with a viewing angle of 45° for the top panel 
and the case of ray 2 with a viewing angle of 0° for the middle 
and bottom panels. The normalisation values are the same as in 
Fig. [131 



The eff'ects of wavelength resolution on the spectral broad- 
ening of the specific intensity profile lx(t) for the base model - 
193 are shown in Fig. [14] As for Fig. [13] here we show the case 
of ray 1 (node) with a viewing angle of 45° for the Doppler ve- 
locity and the case of ray 2 (anti-node) with a viewing angle of 
0° for the line width and the parameter. We fix the spatial 
resolution to OR. Diff'erent spectral resolutions are produced by 
integrating the (high wavelength resolution) intensities over spe- 
cific wavelength bins. We can see that no changes are produced 
in the Doppler velocity. This is due to the lack of asymmetries 
in our model: Regardless of the non-Gaussianity of the profiles 
they will always be symmetric (with respect to the wavelength 
corresponding to the bulk velocity of the plasma). The eff'ects on 
the line width are also minimal, with slightly increasing factors 
between 1.04 and 1.09 among the models. This increase is due 
to the wider wavelength bins over which intensity is added. At 
lower resolution the bumps created by the density profile cannot 
be detected, therefore increasing the Gaussianity of the profile. 
Factors of increase range between 1.20 and 1.23 between the 
models. 



3.4. Differences between models 

The results presented in the previous section apply in a qualita- 
tive manner to all the models presented in this paper. Apart from 
the diff'erences in the integrated intensity between the 171 A line 
(base mod el -1 71) and the 193 A line (base model - 193) treated 
diff'erences between spectral broadening of the 



m section 



3.1 



diff'erent models were also briefly stated in section [X2] and 3.3 



Figure 15 shows the time evolution of the spectral line broad- 
ening (the line width and the non-Gaussianity x^ parameter) re- 
sulting from the Gaussian fit to the specific intensity for all 
the models presented in this study, taking OR as fixed pixel size. 
From left to right the first two panels correspond to the 0° angle 
case, for rays 1 and 2. The last two panels correspond to the 45° 
angle case for rays 1 and 2. The quantities are normalised with 
respect to the temporal average of the corresponding quantity for 
ray 1 with a 0° viewing angle. 
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As can be seen in the panels containing the line profile width 
all models except the long A model have essentially the same 
variation. This is mainly due to the fact that the same velocity 
field is achieved for the same wave number. Longer wavelengths 
(small wave number) produces higher velocity fields, leading 
to larger variation in line widths. Higher velocities also lead to 
stronger compression and rarefaction, and therefore larger vari- 
ation in temperatures and densities. The eff'ects on the spectral 
broadening of the line are low, however, since plasmas with a 
hotter temperature produce lower intensities in the 171 A line. 
The longer wavelengths produce wider regions of compression 
and rarefaction, thus modifying the emission crossed by rays 1 
and 2 between the models. This can be observed in the 45° an- 
gle case for ray 2, leading to stronger variation in line widths 
midx^ values. This eff'ect disappears for a viewing angle of 60°, 
where ray 2 crosses mixed regions. On the other hand, the high 
T model shows in general the lowest values in line width andx^ 
parameter values (within models with same wave number espe- 
cially). This is due to the hot plasma component at rest outside 
the cylinder. Only minimal diff'erences occur between the base 
model - 171 and base model - 193. This is because the emission 
measure weighted radial velocities are very similar between the 
two models. 



4. Discussion 

4.1. Intensity modulation 

Table [2] summarizes the main findings of this work. These re- 
sults have been obtained assuming equilibrium ionisation for the 
considered ions, and should be discussed taking this assumption 
into account. As discussed in section |3.1| based on the results 
in Table [T] in the case of a Maxwellian plasma distribution non- 
equilibrium ionisation severely aff'ects the modulation of inten- 
sity, but does not alter the results obtained for the modulation 
of spectral line shape (width and Gaussianity). The results for 
intensity presented in Table [2] and discussed in this section con- 
stitute therefore upper limits to the more realistic case of plasma 
modulation in a simulation where the ion balance equations are 
solved. 

The modulation of the integrated intensity, especially impor- 
tant for imaging instruments, severely depends on the selected 
spectral line(s) of the observation. In the case of a coronal line 
formed at temperatures below those present in the structure, such 
as the 171 A line for the base model -171, the intensity is only 
lightly modulated by the sausage mode, even when increasing 
the mode's wavelength, with maximum amplitudes of 2 - 4 %. 
Such models present intensity variations in anti-phase with the 
temperature and density variation, due to the negative slope of 
the ion population function at these temperatures. However, the 
modulation can be increased by a factor of 10 - 15 when consid- 
ering coronal lines with formation temperature higher than that 
of the plasma, as is the case for the base model - 79 J. In this case 
we have an in-phase variation between these quantities. 

In our model the intensity is governed by the ion popula- 
tion fraction rather than the density of emitting ions. Indeed, as 
shown by Figure |6j the change in the contribution function from 
the temperature modulation (at average density) is one order of 
magnitude larger than the change from the density modulation 
(at average temperature). This implies that non-equilibrium ioni- 
sation eff'ects will be important, especially for Fe xii, in which the 
modulation of intensity would completely disappear for the short 
timescales of the sausage mode (as discussed in section[3]T]). The 
important role of the contribution functions of spectral lines. 



especially concerning non-equilibrium ionisation eff'ects, have 
been discussed by |De Moortel & B radshaw (|2008j). Through 
synthetic intensity modelling of TRACE channels, the authors 
show similar results as those stated above. They further show 
that double periodicities can be obtained from the shape of the 
response function of the instrument due to perturbations leading 
to temperatures that drift in and out of the maximum formation 
temperature of the line. 

Variation of the modulation of the intensity by kink and 
sausage modes due to variation of the angle between the line- 
of-sight and the emitting struc ture were first analysed analyti- 
cally by |Cooper et al.| ( [2Q03a|b| ). Considering the eff'ect of curva- 
ture they showed that propagating kink and sausage modes could 
produce an increas e of intensity along part o f the loop, such as 
that observed by Willi ams et"aL] ( |200 1 1 [2002| ) . In their work the 
intensities are approximated by the line-of- sight column t hick- 
ness / and the density p, as / p^l. As discussed in section 3.1 



the neglect of the contribution function can be misleading, anc 
the variation of line intensities can be severely underestimated. 
In the present work we have monotonically decreasing depen- 
dence on the line-of- sight angle for the range of parameters con- 
sidered here, match ing the results of [Cooper et al.| ( [2003a l. In 
Cooper et al. ([2003b ) the case of a smooth density profile is con 



sidered (with the corresponding dispersion relations associated 
with that profile), which introduces an optimal angle for which 
the intensity variations are largest. This is due to the fact that in a 
2D model the MHD modes are more compressive, and therefore 
produce significant density variations outside the cylinder. This 
eff'ect is however severely reduced for 3D geometries, as shown 
by |Luna et al. ( 2008| ), suggesting that the existence of an optimal 
angle could be debated. The eff'ect of such density profiles in a 
3D model such as the one considered here is the subject of future 
research. 

[ Gruszecki et al.| p012| ) have performed 3D numerical sim- 
ulations of (trapped and leaky) sausage modes oscillating in a 
cylinder with a smooth density profile and a nalysed the modu - 
lation on the integrated intensity. As in Cooper et al.| ( |2003a[ ), 
the intensity is approximated by integrating the density squared 
along the line-of-sight (for a viewing angle of 0°), but they fur- 
ther consider the pixel size resolution. With the same amplitude 
as used in this work they find intensity variations on the order of 
2 % with OR pixel size resolution, which are strongly reduced by 
a factor of roughly 15 when increasing the pixel size to 3R (one 
wavelength), thus matching the results presented here. As can 
be seen in the table, among all parameters, we find the pixel size 
resolution to have a dominant role. A pixel size resolution on the 
order of the wavelength erases the longitudinal structuring in all 
models. 

The sausage mode in the corona oscillates with relatively 
short periods. Indeed, in the non-leaky regime and for low-yS 
plasma conditions, the period of an MHD sausage mode can 
be shown to satisfy the relation P < where 70 ~ 2.4 



is the first zero of the Bessel function J{){x) ( [Nakariakov et al.| 
[2003 ; Ko pylova et al . |2007 ) . In our models, the periods oscillate 
around 5 s, making the cadence of the instrument a very impor- 
tant parameter. Intensity variations can be decreased by a factor 
of 1.5 and more than 3 to 5 periods of the mode need to be ob- 
served in order to half resolve the full profile. 

All the interpretations as sausage modes in observations are 
mostly based on qualitative arguments, especially restricted by 
the spatial resolutions of the instrument, which in most cases 
does not allow to investigate in detail the longitudinal varia- 
tion of the oscillations along the flaring coronal structure. Still, 
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Table 2. Influence of line-of-sight angle, and spatial, temporal and wavelength resolution on the modulation of observable quantities 
caused by the sausage mode. 



Quantity 


Ray 


Amplitude 


Angle 

0° 60° 


Spatial 
R^ 3 R 


Cadence 
s^ 2 s 


/l(mA) 
1.4 8.5 


Period 


Intensity 


1 


~0 


~0 


~0 


i 1.2-1.5 




P 


2 


2 - 4 % (35 %, 193) 


i2.5 


UO- 17 






Doppler 


1 


< 1 km s"^ 


1 


1 


i 1.1 - 1.4 


1 


P 


velocity 


2 


7 km s"^ 


I* 


1 




Line 


1 


7 - 20 % 


I* 


T* 


i 4.6 - 6.3 


i 1.04-1.09 


p 


width 


2 


38 - 49 % 


i 5.6 -11.3 


i 2.5 -3.4 


2 




1 


2.1-2.9 




I * 


i4-5 


T 1.20-1.23 


P P 


2 


3.75-4 


i 2.1 -2.5 


12.2-5.5 


4' 2 



Notes. Since the sausage mode introduces longitudinal structuring we show values for rays 1 (node) and 2 (anti-node) for geometrical quantities. 
The amplitude denotes the maximum variation in time with respect to the value at t=0 of the quantity considering all parameters. The unit 
corresponds to percentage except for the Doppler velocity where it denotes the maximum value, and the parameter, where it denotes the 
variation in orders of magnitude. A range in the values of the quantity corresponds to variations between the models. The values in the other 
columns denote the factor by which the variation of the quantity is modified when varying the specific parameter of the column. The direction 
of the parameter variation is indicated in the column head. Upside and downside arrows (or both) denote, respectively, an increase or a decrease 
(or both if non-monotonic) by the specified factor as the parameter is varied. An asterisk indicates that the variation in time of the quantity is 
negligible for the initial value of the parameter. Values in parenthesis correspond to some relevant values for the base model - 193. The last column 
corresponds to the possible observed period in the specific observed quantity, with P the period of the sausage mode. 



resolutions such as that of NoRH have allowed the detection 
of difl'erent oscillatory power between loop top and footpoint, 
and allow more precise identification of the mode harmonics 



(INakariakov et al.|2003||Melnikov et al.|2005l|Inglis et al.|2008| 



[Inglis & Na kariakov[2Q09). In these works, the fundamental 
mode and sometimes the first overtone can be identified. As dis- 



cussed in Andries et al. ( 2009 ), observations of higher overtones 
are specifically important, since they significantly improve phys- 
ical parameter estimates from application of MHD seismology. 
Observations of a fundamental mode at the apex of loops would 
be comparable to the case of ray 2 in our model, while observa- 
tions towards the footpoints would correspond better to the case 
of ray 1 . Most observations of sausage modes report short and 
thick loops, with lengths around 25 Mm. In this case, the res- 
olution of NoRH (5'' - W) corresponds to a pixel size below 
y or 1^, with ^0 and being the wavelengths of the funda- 
mental mode and the first overtone. According to the results pre- 
sented here a pixel size of |^ reduces the intensity variations by 
a factor of 30 % - 40 %. For detection of higher harmonics the 
pixel size starts to severely decrease the intensity variation, as 
seen in the lower bottom panel of Fig. [5] Observations of modes 
with large line-of-sight angles with respect to the normal to the 
loop plane can easily have rays crossing mixed regions of expan- 
sion and contraction, corresponding for instance to nodes of the 
first overtone and the fundamental mode. These scenarios would 
bring results that are qualitatively similar as those obtained in 
the present model when varying the line-of-sight angle. 



4.2. Spectral line broadening modulation 

Table [2] especially shows that even when intensity variations are 
low (as in the base model -171), spectral line broadening can be 
important, and could therefore be observed by current spectrom- 
eters. In section 32^ we have shown that the nature of the spectral 
line broadening from the sausage mode is not thermal but mainly 
turbulent broadening. This implies that non-equilibrium ionisa- 
tion will not influence the obtained spectrometric results. 

The sausage mode is essentially a radial mode, and therefore 
Doppler motions (through single-Gaussian fitting) c an o nly be 
detected when looking at an angle as shown in sectionlT2] In this 



case, bulk motions can be detected as blueshift and redshift ex- 
cursions below 7 km s"\ roughly one third of the maximum ve- 
locities achieved in the present model. Maximum variations are 
obtained for 0° viewing angle for the line width and the pa- 
rameter for ray 2 (anti-node), with values close to 50 % and 4 or- 
ders of magnitude, respectively (see table). As shown in Fig.pT 
longitudinal structuring can be detected in Doppler velocity, line 
width and the;^^ parameter for angles up to 30°, pixel size res- 
olutions up to IR, temporal resolution of 1 s (although with 
2 s resolution it is still possible to detect bulk Doppler veloci- 
ties) and for all considered wavelength resolutions. Furthermore, 
the line width and x^ parameter present double periodicities. 
For lower resolutions or higher viewing angle the parameter 
presents quadruple periodicity, but in these cases all quantities 
exhibit variations decreased by factors of 2 to 10 for each pa- 
rameter, and therefore may become barely detectable. 

So far very few reports exist of periodic spectral broaden- 
ings in EUV lines with little or no associated intensity varia- 
tions interpreted as fast sausage modes. Most of these st udies are 
performe d with the Hinode/EIS imaging spectrometer ((Culhane 
etaT72007), with the slit (with a pixel size resolution of 2'') 
and tempora l resolutions on the order of 10 s or more. Kita gawa| 
et al.| ( 2QTQ| Fourier analyze EUV intensity, Doppler velocity 
and line width time series of data in an active region observed in 
the Fexii 195 line. Moss regions were found to exhibit signifi- 
cant correlated intensity and Doppler velocity oscillations, with 
variations of 0.4 % - 5.7 % and 0.2 - 1.2 km s"^ respectively, 
consistent with fast MHD waves modulation. Depending on the 
significance level of intensity oscillations to Doppler velocity 
and line width oscillations the waves were interpreted as either 
fast sausage modes, Alfven or kink modes. However, our results 
suggest that a low significant level in intensity can be obtained 
with a sausage mode (any pixel size resolution for 171 line, 3R 
pixel size resolution for 193), while keeping a relatively high sig- 
nificant level for Doppler velocities and line width (above 10 % 
variation). The long reported periods on the order of minutes 
could be produced by leaky fast sausage modes, especially to- 
wards the footpoints of loops, as explained below. 

Many reports exist of observations in the same lines hav- 
ing similar attributes, interpreted as slow modes, kink modes 
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De Moortel & Nakariakov 



or Alfven waves (see review by 
2Q12| ). Although the authors discard the sausage mode interpre- 
tation for other reasons, the essential characteristic discarding 
sausage modes are the long reported periods. Sausage modes 
with long periods can exist in the corona, as shown by Pascoe, 
et al.| ( [2n07b| ), but are leaky and therefore their observability is 



severely limited for ranges matching the reported periods. 

Often, observations of propagating disturbances are inter- 
preted as slow modes based on the phase speed closely match- 
ing the local sound speed, periods between 2 and 10 minutes, 
(low) intensity variations (of a few percent with respect to the 
background intensity) and the intensity fluctuation lagging the 
Doppler shifts by 1/4 period ( jWang et al.|2Q 03). Here we would 
like to stress the fact that it is the combination of all these char- 
acteristics that supports such interpretation. Apart from the first 
two, (phase speed and period), a sausage mode can also explain 
the rest (intensity modulation and phase shift, as obtained when 
observing with a non 0° angle). 

Wavelength resolution is found to play a minor role, relative 
to spatial and temporal resolutions. As shown in section [34] only 
minor diff'erences appear in the spectral broadening (and none 
in the Doppler velocity) when decreasing the wavelength reso- 
lution to only 8 points across the spectrum (AA =8.5 mA). This 
may be due to the symmetric nature of the sausage mode (with 
respect to the axis of the cylinder), but also to the absence of 
non-uniformities in density or temperature in the present mod- 
els. Pascoe et al. ( 2QQ7a| 2QQ9| ) considered the eff'ects of varying 
loop cross-section and fine multilayered structuring on the reso- 
nant periods and the spatial profiles of fast sausage mode oscil- 
lations. These numerical studies find that the resonant properties 
of long-wavelength sausage standing modes of the slab are not 
sensitive to the details of fine structuring and the spatial profiles 
are only weakly aff'ected by the cross- section inhomogeneity. 
This suggests that the inherent symmetries of the fast sausage 
mode are rather robust to loop divergence and internal density 
fine-structuring, and therefore that the properties of the modula- 
tion found for simple cylindrical geometries, such as those used 
here, are robust as well. 

As discussed in the previous section, the detection of sausage 
modes is severely limited by its leaky nature in most coronal 
conditions. On the other hand the lack of reported observations 
of sausage modes with spectrometers may also be due to an in- 
strumental reason. Our results show that high spatial and tempo- 
ral resolution in EUV lines is essential for the detection of such 
waves. Perhaps the most promising instrument for the detection 
of these waves in the near future is IRIS, the Interface Region 
Imaging Spectrograph, whose launch is scheduled for the spring 
2013. IRIS is specially designed for spectrometric observations 
of the solar atmosphere in EUV lines with high spatial and tem- 
poral resolutions. For instance, observation in the Fexii 1349 A 
line (log r = 6.2 K) will have a cadence of 1 - 2 s, a spatial res- 
olution of 0.4'' (with a field of view of 0.3'' x 40"), an eff'ective 
area of 2.8 cm^ and a wavelength resolution of AA = 12.5 mA. 

The new Atacama Large Millimeter/submillimeter Array 
(ALMA) telescope has started to operate and off'ers unique ad- 
vantages for solar observations. Spatial resolutions from 0.06" 
at 85 GHz to 0.005" at 950 GHz and short integration times of 
32 ms are expected, although for solar observations a special fil- 
ter is needed, which may considerably reduce the capabilities 
of the telescope. Nevertheless, flare and chromospheric research 
may considerably benefit from ALMA, especially in the subject 
of fast waves such as the ones studied here. 



The results obtained here would correspond to the ideal case 
of an imaging spectrometer with at least 2 s cadence, a pixel size 
resolution of at least one wavelength, and a wavelength resolu- 
tion of at least 8 points across the line. Furthermore no blending 
of lines emission was considered, as is produced when observ- 
ing with the channels of the AIA instrument on board the Solar 
Dynamics Observatory, no background optically thick contribu- 
tion and no instrumental noise. Together with the eff'ects of cur- 
vature and non uniform media (leading to diff'erent wave disper- 
sion relations) these scenarios will be treated in future work. 

5. Conclusions 

In this paper we aimed to understand the variations of observ- 
able quantities introduced by the fast MHD sausage mode in 
coronal structures, considering geometrical and instrumental ef- 
fects. Complexities introduced by non-uniformity of the cylinder 
or the background, curved geometries, or complex instrumen- 
tal eff'ects such as instrumental noise are not taken into account. 
Although important for realistic models of loops, such charac- 
teristics are secondary and their eff'ects on intensity modulation 
cannot be properly tackled before understanding the eff'ects of 
more simple parameters on simple geometrical models such as 
the ones considered in this work. 

Details from the main results of this work can be found in 
Table [2] The main findings can be summarized as follows. 

The intensity modulation is strongly dependent on the forma- 
tion temperature of the spectral line with respect to the temper- 
ature of the plasma. Assuming equilibrium ionisation, for plas- 
mas with higher temperature than the formation temperature of 
the line, such as the case of the base model -171, the modu- 
lation is low (2 % - 4 %) and is slightly aff'ected by the geo- 
metrical variation of the cylinder. On the other hand, for plas- 
mas with lower temperature than the formation temperature of 
the line, such as the base model -193, the modulation matches 
the thermodynamical evolution (density and temperature), and 
can therefore be significant (35 %). Non-equilibrium ionisation 
would completely eliminate the Fe xii intensity modulation (and 
reduce by a factor of 1.5 or more the Feix intensity) due to the 
fast timescales of the sausage mode. Simple approximations to 
the intensity taking into account the column thickness and dis- 
regarding the contribution function G^QiT^Ue) can therefore be 
very misleading. 

Depending on the contribution function, the intensity can be 
in anti-phase with the density (or temperature) profile (case of 
the base model - 171) or in-phase (case of the base model - 193). 
The line width and the;^^ parameter are in phase with the square 
of the radial velocity. The parameter can also respond to the 
density variation, which becomes significant for coarse pixel size 
or large viewing angle. The Doppler velocity is in anti-phase 
with the radial velocity when watching at an angle (non 0°). 
Correspondingly, for these angles, the Doppler velocity presents 
a phase shift of one fourth of the period with the intensity. In 
accordance with the previous result, the intensity and Doppler 
velocity exhibit the same periodicity as the mode. Double peri- 
odicities are characteristic of the line width and the;^^ parameter. 
The parameter can present quadruple periodicities for coarse 
pixel size or large viewing angle. 

The detection of intensity modulation from sausage modes is 
strongly dependent on the angle of observation. The intensity de- 
creases with the viewing angle (where 0° angle denotes the per- 
pendicular to the axis of the cylinder) by a factor of 2.5 for the 
studied range. Furthermore, the intensity modulation can vary 
significantly longitudinally (along the axis of the cylinder). Rays 
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crossing the cylinder at nodes exhibit considerably less variation 
than those crossing the anti-nodes (a factor of 20 or more, de- 
pending on the pixel size resolution). 

Setting aside equilibrium ionisation effects, the detection of 
intensity modulation is strongly dependent on spatial and tempo- 
ral resolutions. Both can be the dominant parameters, especially 
if the pixel size is on the order of the mode's wavelength or the 
temporal resolution is on the order of the mode's period. A pixel 
size or a cadence on these orders (respectively) can completely 
erase the longitudinal structuring, therefore eliminating the pos- 
sibility of detection. Wavelength resolution is not particularly 
important for the case of sausage modulation. This is due to the 
inherent symmetry of the sausage mode with respect to the axis 
of the cylinder, leading to symmetric line profiles. Numerical 
studies suggest that this as well as most results presented here 
may be robust to the presence of inhomogeneities in the model 
(leading to possible assymmetries), as well as effects from vary- 
ing cross- section along the cylinder. 

Even when the intensity modulation is low, the variation of 
the Doppler velocity, and especially the line width and the 
non-Gaussianity parameter can be significant (a factor of up to 
40 % for line width and up to four orders of magnitude change 
for the;^^ parameter), placing spectrometers in a clear advantage 
over imaging instruments for the detection of these waves. This 
result is independent of the ionisation state of the plasma since 
the spectral line is mainly affected by the unresolved turbulent 
motions rather than the thermal motions. Significant variability 
in all quantities can still be obtained when viewing at an angle 
of up to 30° (Doppler velocities can especially be detected for 
the range 15° - 45°), with pixel size resolutions up to IR, or 
temporal resolution of 1 s (although with 2 s resolution it is still 
possible to detect significant Doppler velocities) and for all con- 
sidered wavelength resolutions. We consider this to be the most 
important result of this paper. 
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